pub fn new_birthday_probability(n: u32) -> f64 {
    1f64 - (1..n + 1).map(|n| (365 - n + 1) as f64 / 365f64).product::<f64>()
}
